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On the basis of the multi-orbital dynamical mean field theory, a three-orbital Hubbard model 
with a relativistic spin-orbit coupling (SOC) is studied at five electrons per site. The numerical 
calculations are performed by employing the continuous-time quantum Monte Carlo (CTQMC) 
method based on the strong coupling expansion. We find that appropriately choosing bases, i.e., 
the maximally spin-orbit-entangled bases, drastically improve the sign problem in the CTQMC 
calculations, which enables us to treat exactly the full Hund’s coupling and pair hopping terms. 
This improvement is also essential to reach at low temperatures for a large SOC region where the 
SOC most significantly affects the electronic structure. We show that a metal-insulator transition 
is induced by the SOC for fixed Coulomb interactions. The insulating state for smaller Coulomb 
interactions is antiferromagnetically ordered with the local effective total angular momentum j = 
1/2, in which the j = 1/2 based band is essentially half-filled while the j = 3/2 based bands are 
completely occupied. More interestingly, for larger Coulomb interactions, we find that an excitonic 
insulating state emerges, where the condensation of an electron-hole pair in the j = 1/2 and j — 3/2 
based bands occurs. The origin of the excitonic insulator as well as the experimental implication is 
discussed. 

PACS numbers: 71.27.+a, 71.30.+h, 75.25.Dk 


Recent experiments have reported interesting obser¬ 
vations for 5 d transition metal Ir oxides in the layered 
perovskite structure such as S^IrCU and Ba 2 lr 04 [1-5]. 
In these materials, along with moderate electron corre¬ 
lations [6, 7], there exists a strong relativistic spin-orbit 
coupling (SOC), which splits t^ g orbitals, already sepa¬ 
rated from e g orbitals due to a large crystal field, into 
the effective total angular momentum j = 1/2 doublet 
and j = 3/2 quartet orbitals in the atomic limit [8]. 
Since there are nominally five 5 d electrons per Ir ion, the 
j = 1/2 orbital is half filled while the j = 3/2 orbitals 
are fully occupied. 

As opposed to simple expectation from strongly cor¬ 
related 3d and 4 d transition metal oxides [9-13], the ex¬ 
periments have revealed that the ground state of these Ir 
oxides is a j = 1/2 antiferromagnetic (AF) insulator [14- 
17]. The theoretical understanding of the j = 1/2 AF 
insulator has been also reported [14, 18-24]. Moreover, 
even possible unconventional superconductivity has been 
proposed once mobile carriers are introduced into the in¬ 
sulating state [25-28]. However, the electronic structure 
in multi-orbital systems with the competition between 
the electron correlations and the SOC has not been thor¬ 
oughly understood. When the SOC is significantly large, 
the j = 1/2 based band is completely separated from 
the j = 3/2 based bands, and thus a single-orbital de¬ 
scription of the j = 1/2 based band is expected to be 
valid. On the other hand, when the SOC is small, this 
picture breaks down and the electronic structure should 
be largely affected not only by Coulomb interactions but 


also by the multi-orbital nature. Therefore, a question 
naturally arises: what is the ground state of multi-orbital 
systems with the competition between the electron cor¬ 
relations and the SOC. This is precisely the main issue 
of this paper and we will demonstrate the emergence of 
exotic insulators in a three-orbital Hubbard model with 
the SOC. 

Here, on the basis of the multi-orbital dynamical mean 
field theory (DMFT) [29], we numerically study a three- 
orbital Hubbard model with the SOC at five electrons 
per site, corresponding to ( t 2 g ) 5 electronic configuration. 
The continuous-time quantum Monte Carlo (CTQMC) 
method based on the strong coupling expansion is em¬ 
ployed as a multi-orbital impurity solver [30]. We find 
that the sign problem is significantly improved by ap¬ 
propriately choosing the maximally spin-orbit-entangled 
bases, which allows us to treat exactly the full Hund’s 
coupling and pair hopping terms. This improvement is 
crucial also for the calculations in a low temperature and 
a strong SOC regions where the metal-insulator transi¬ 
tion (MIT) occurs. We find that for fixed Coulomb in¬ 
teractions the MIT is induced by increasing the SOC. 
The insulating phases include the j = 1/2 antiferromag¬ 
netically ordered phase as well as a multi-orbital AF in¬ 
sulating (AFI) phase. In addition, we find an excitonic 
insulating (EXI) phase, where an electron-hole pair in 
the j = 1/2 and j = 3/2 based bands are condensed. 

The three-orbital Hubbard model studied 
here is described by the following Hamiltonian: 
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H = H 0 + Hi, where H 0 = £ (M ,> £ 7>CT ^4^'^ ~ 

vY,i,-y,a n lcT + AE il 7 ,w<7|L^> • (o-|Si|f 7 ')ct 7 (T c i5(T / 
represents the non-interacting part of the model 

and Hi = C/E i) 7 «7t n 4 + Ul f 1 T,i n ^a n la n t + 

J2i, 17 tS,a n l<r n ia ~ J Ei, 77 t5 + 

J' J2i, j7 is c l 7 1 c 1 7 4 c i< 54 c i <51 describes the local Coulomb 
interactions. Here, t 1 sets the nearest-neighbor hopping 
amplitude for t^g orbitals 7 = (d yz ,d zx ,d xy ) on the 
Bethe lattice with coordination number Z [31]. We 
adopt a semielliptic density of states for Z 00 
assuming the same bandwidth for the three orbitals 
(t 7 = t/VZ), where the DMFT is exact [31, 32]. 
c\ ( Ci 1(7 ) is an electron creation (annihilation) op¬ 

erator with spin a (=t, i) and orbital 7 at site i and 
n la = C Jjo- c i 7 < 7 - ^ is SOC and Li (Si) is the orbital 
(spin) angular momentum operator at site i. The 
chemical potential /i is tuned to be at five electrons per 
site, a denotes the opposite spin of a. Hi includes the 
intra (inter) orbital Coulomb interaction U ( U '), the 
Hund’s coupling (J), and the pair hopping (J'). We set 
U = U' + 2 J and J = J' = 0.15 U [33]. Employing the 
CTQMC method as an impurity solver [30] to calculate 
the imaginary-time Green’s functions at the impurity 
site, G^(i,r) = -(TVci 70 -('r)c| 5 <T ,( 0 )), we can solve 
numerically exactly the model. In what follows, C7, 
A, temperature T, and frequency uo are in units of t. 
We also omit the site index in the Green’s function, 
(t) = (i, r), unless G^ (i, r) is site dependent. 

Let us first examine the numerical accuracy of the 
CTQMC calculations. Generally, a negative sign prob¬ 
lem is one of the most serious issues in QMC calcula¬ 
tions. Specially, when QMC methods are employed for 
the multi-orbital DMFT calculations, the sign problem 
seriously prevents us from simulating at low tempera¬ 
tures, which often enforces to approximate the Hund’s 
coupling and pair hopping terms [34]. Indeed, for our 
model, we find that the sign problem becomes destruc¬ 
tively serious at low temperatures specially when the 
SOC is large (see Fig. 1 ). 

It is now important to recall that the sign problem 



FIG. 1: (Color online) (a) A dependence of the average sign 
(sgn) for U = 8 at T — 0.08. (b) T dependence of (sgn) for 
U = 8 at A = 0.24. For both t 2 g and j bases in the CTQMC 
calculations, the paramagnetic and orbital disordered solu¬ 
tions are assumed. 


in QMC calculations is a basis dependent problem and 
it can be improved by appropriately choosing bases [35]. 
Although we encounter the serious sign problem when the 
original t^g bases (ci la ) are used, we find that the max¬ 
imally spin-orbit-entangled j bases (a^- m ) improve the 
sign problem significantly. These j bases are the eigen¬ 
states of Hq in the atomic limit and are related to the 
t 2 g bases, 

>/2 -iy/2s y/2s \ / c idyzd - \ 

£ i 2 J I Cid zx Q- J ,(1) 

—>/3s —W 3 0 / \ Cid xy a J 

where s = 1(—1) for a =t (i). In the j bases represen¬ 
tation, the number of off-diagonal elements in G^° a (r) is 
reduced and as a result the sign problem can be allevi¬ 
ated. 

Figure 1 shows typical results of the average sign (sgn) 
(see Ref. 35 for the definition) for V = 8 calculated for 
the paramagnetic and orbital disordered solutions. The 
A dependence of (sgn) at T = 0.08 in Fig. 1(a) clearly 
shows that (sgn) ~ 1 for the j bases, much larger than 
(sgn) for the t<i g bases particularly for large A. Simi¬ 
larly, the T dependence of (sgn) at A = 0.24 in Fig. 1(b) 
demonstrates the remarkable improvement of (sgn) for 
low T when the j bases are used. These improvements 
of (sgn) guarantee higher accuracy of our CTQMC cal¬ 
culations with less computational cost in a wide range of 
A even at very low T. 

First, we shall briefly examine the MIT induced by the 
SOC in paramagnetic and orbital disordered states. A 
typical example of the evolution of single-particle excita¬ 
tion spectrum Aj ?m (cj) = — 1 /MmG^ (uj + i0 + ) in the 
real frequency uo with increasing A is shown in Fig. 2(a). 
Here, 0 + is positive infinitesimal. The maximum en¬ 
tropy method is employed to calculate Aj^ m (uj) from 
Gj’JJ(r) [36]. As shown in Fig. 2 (a) for U = 8 , the 
SOC induces the transition from the metallic to insulat¬ 
ing state at A ^ 0.31, where the quasi-particle peak near 
the Fermi level vanishes and the single-particle excitation 
gap starts to open. The insulating state for A > 0.31 is 
the j = 1/2 Mott insulator, where the j = 1/2 based 
band is half-filled while the j = 3/2 based bands are 
fully occupied [see Fig. 2(b)]. Figure 2(b) also indicates 
that the j = 1/2 based band, which is degenerate with 
the j = 3/2 based bands at A = 0 , gradually separates 
form the j = 3/2 based bands with increasing A [see also 
Fig. 2 (c)]. This is essential for the SOC induced Mott 
transition because the critical U for the MIT decreases 
as the orbital degeneracy is degraded [37]. 

Now, we shall consider possible ordered states. To this 
end, we introduce magnetic order parameters, Mj jTn (l) = 
^T l m'=±m si 9 n(m')(al jm ,ai jm >), where l(= A,B ) indi- 
cates two sublattices [38] and aij m is defined in Eq. (1) 
with site i being on sublattice /. In addition, we in¬ 
vestigate excitonic orders formed by an electron-hole 
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FIG. 2: (Color online) (a) Single-particle excitation 

spectrum Aj jTn (uj) and (b) electron density nj jTn = 
Y2 m '=±m( a lj m ' a ijm') for U = 8 at T = 0.06 with varying 
A. Paramagnetic and orbital disorder states are assumed, (c) 
Schematic density of states for the non-interacting limit with 
A = 0 and A ^ 0. Fermi level is at uj = 0. The j = 3/2 based 
bands are completely filled when A > 4/3. 


pair in different (j, m) bands with excitonic order pa- 
rameters, ^ (l) = where (j,m) ^ 

(j',m'). We also calculate the electron density rij ?m (/) = 
J2m'=±m( a ljm' a i-jrn'), Aj >m (l,u), and the anomalous 
excitation spectrum F 3 ^ (/, uj)=— l/nlrnGj (1 ,uj + 
i0 + ) [39, 40]. Since we have found that these order pa¬ 
rameters satisfy Mj jTn (A ) = — My m (T>) and /S 3 - (A) = 

A i-m ( B )i as wel1 as n j,rn( A ) = n j,m{ B )i A j,m { A = 
Aj-miB.u) and F 3 ^ 1 (A,c j) = F 3 F^{B,uj), we will 
drop the sublattice index l in these quantities hereafter. 

Let us first explore a case with U = 8. Figure 3 (a) 
shows the A dependence of Mj ?m , Aj , and n J?m . First 
of all, it is apparent that there is no excitonic order for 
all values of A. Second, for A > 0.25, only Mi/ 2 , 1/2 is 
finite with n 1/2 , 1/2 ~ 1 and n 3 / 2 , 1/2 = n 3 / 2 , 3/2 ~ 2, 
revealing the j = 1/2 AFI state. Third, for A < 0.25, the 
magnetic order disappears and n 3 / 2?1 / 2 = n 3 / 2 3 / 2 starts 
to decrease from two with decreasing A, implying the 
breakdown of the single-orbital description of the half- 
filled j = 1/2 based band. 

Figure 3(b) shows Aj iTn {uj) for two A’s at U = 8. Sim¬ 
ilar to the paramagnetic cases (see Fig. 2), for large 
A = 0.32 in the AFI phase, a finite gap is clearly open 
in the j = 1/2 based band while the j = 3/2 based 
bands are fully occupied (i.e., band insulators). How¬ 
ever, notice that a i/ 2 ,i/ 2 (u) [ a i/ 2 ,-i/ 2 (u)] has larger 
(smaller) weight in the occupied states than in the un¬ 
occupied states for sublattice A because of the AF or¬ 
der. For smaller A = 0.08, besides broad structures 
around uj ^ ±U/2 , the sharp quasi-particle peak appears 
around the Fermi level for both j based bands, indicating 
a metallic state with strong correlations. 

Next, we shall examine a case with larger U = 9.25. 
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FIG. 3: (Color online) (a) and (c): A dependence of stag¬ 
gered magnetization Mj,m (top), excitonic order parameter 
A J j’ ™ (center), and electron density (bottom), (b) and 
(d): Single-particle excitations spectrum Aj,m(uj) for two A’s. 
Other parameters used are T = 0.06 and U — 8 (a,b) and 
9.25 (c,d). Anomalous excitation spectrum ^f /2 [± 1/2 (aO is 
also shown in (d). 


As shown in Fig. 3(c), the j = 1/2 AFI phase is found 
with no exitonic order for A > 0.25, the same phase dis¬ 
cussed above for smaller U. Aj :7n (uj ) shown in Fig. 3(d) 
clearly exhibits a finite gap only in the j = 1/2 based 
band while the j = 3/2 based bands are fully occupied. 
The difference appears as compared with the case with 
smaller U when A decreases. For 0.12 < A <0.25, all 
three Mj ?m ’s are now finite and ni/ 2j i/ 2 (^ 3 / 2 , 1/2 and 
713 / 2 , 3 / 2 ) increases (decrease) from one (two) with de¬ 
creasing A [seen Fig. 3(c)]. We further confirm a finite gap 
in Aj^ m (uj) for both j based bands. These results indicate 
that this phase is a AFI state but apparently breaks the 
single-orbital description of the half-filled j = 1/2 based 
band. 


More remarkably, for A < 0.12, we find that the exci¬ 
tonic order parameter A^’^i/^ is finite with non zero 
AF order [see Fig. 3(c)]. Although the magnetic order 
eventually disappears for A < 0.04, A^’ii /2 remains fi¬ 


nite. The most important feature in this phase is that 

A 3/2,1/2 . A 3/2,-1/2 _ n , , A 3/2,1/2 A 3/2,-1/2 . ~ 
^1/2,1/2 ^1/2, — 1/2 — 0 but A-j_/ 2jl / 2 ; ^ 1 / 2 , -1/2 1S ^ 


nite [41] and staggered between the two sublattices, indi¬ 
cating that the excitonic order is accompanied by trans¬ 
lational symmetry breaking. It is also noticed in Fig. 3(c) 
that n 3 / 2?1 / 2 7 ^ n 3 / 2?3 / 2 in this phase because of the pres¬ 
ence of non-zero excitonic order. Moreover, we find that 


a single-particle excitation gap is always finite in both 






















































4 


(a) 



(b) 

0.4 

0.3 

^ 0.2 
0.1 
0 


j=1/2 AFI 


MOAFI 


Metal E^I+AF^ 

; EXI + 

9 001 10 

U 


FIG. 4: (Color online) (a) U dependence of the effective SOC 
A* (see the text for definition) for A = 0.02 at T — 0.06. For 
comparison, the results for A = 0.08 and 0.16 are also shown, 
(b) U -A phase diagram at T — 0.06. (MO)AFI, EXI, and 
OOI stand for (multi-orbital) antiferromagnetic insulating, 
excitonic insulating, and orbital ordered insulating phases, 
respectively. A plus mark indicates a set of U and A where 
the effect of a tetragonal crystal held is examined. 


j based bands, including the phase where only the exci¬ 
tonic order parameter is non zero, as shown in Fig. 3(d). 
In addition, in the EXI phase, ^f/ 2 ’± 1/2 ( CJ ) are f° un d 
to be non zero. This is a strong evidence for the exis¬ 
tence of a stable EXI state formed by an electron-hole 
pair between the j = 1/2 and j =3/2 based bands with 
m = ±1/2. In the limit of A = 0, this EXI state is re¬ 
placed by the orbital ordered insulating state in the t^g 
orbitals reported previously [42]. 

We shall now argue that there are two important in¬ 
gredients for the emergence of the EXI state. First, A is 
strongly renormalized by the Coulomb interactions and 
the renormalized A becomes sizably large even for small 
A in the metallic phase at the vicinity of the EXI phase. 
To show this, we calculate the first moment of A/ ?m ( cj), 
defined as Wj :Tn = f dcu and evaluate the ef¬ 

fective SOC A* = Wi/2,1/2 — Ws/2,1/2 [43]. The U depen¬ 
dence of A* for the metallic phase is shown in Fig. 4(a). 
Indeed, A*/A is strongly renormalized with increasing U 
and can be significantly large, as large as ^ 9 at V = 8.75 
near the transition to the excitonic insulator. 

The other key factors are the Hund’s coupling J and 
the pair-hopping J'. We have performed the static mean- 
filed analysis and found that the EXI state emerges only 
for a small A region with non zero magnetic order. We 
have also found that the excitonic order is never stabi¬ 
lized when J = J' — 0. This implies that the particle 
number fluctuations are essential for the EXI state since 
the particle number of each j = 1/2 and j = 3/2 based 
bands is conserved separately when J = J' = 0. These 
results suggest that not only the SOC but also the full 
Hund’s coupling and the pair hopping play the important 
role for stabilizing the EXI state. The similar tendency 
is found in the DMFT calculations. 

Finally, we have also examined the stability of the EXI 


state [see Fig. 4(b)] in the case where the degeneracy of 
t 2 g orbitals is lifted due to a tetragonal crystal field, caus¬ 
ing a finite splitting S between the d xy orbital and the 
other two orbitals. We have considered both magnetic 
and excitonic orders in the DMFT calculations [44]. Al¬ 
though the crystal field splitting reduces the exciton or¬ 
der parameter as compared to the one at S = 0 , we have 
found that the EXI state is robust against the crystal field 
splitting as long as \S\ is smaller than the single-particle 
excitation gap at S = 0 . 

In summary, we have studied the three-orbital Hub¬ 
bard model with the SOC using the multi-orbital DMFT. 
We have employed the CTQMC method based on the 
strong coupling expansion and found that the sign prob¬ 
lem is significantly improved by using the maximally 
spin-orbit-entangled j bases. The improvement is essen¬ 
tial to treat exactly the full Hund’s coupling and pair 
hopping terms and also to reach at low enough tempera¬ 
tures for large SOC. 

We have applied this method to determine the finite 
temperature phase diagram on the V -A plane at five elec¬ 
tron per site. The U -A phase diagram at the lowest tem¬ 
perature T = 0.06 is summarized in Fig. 4(b). For small 
[/(< 8 ), we have demonstrated that the SOC induces the 
transition from the metallic to the j = 1/2 AFI state, in 
good agreement with the previous numerical study [ 20 ]. 
For large [/(> 8 ), the multi-orbital AFI phase appears for 
intermediate values of A, where the single-orbital descrip¬ 
tion breaks down. Moreover, we have found that further 
overlap of the non-interacting j = 1/2 and j — 3/2 based 
bands favors the EXI phase, either with or without mag¬ 
netic order, which dominates the phase diagram for small 
A (< 0.1). The EXI state is formed by an electron-hole 
pair in the j = 1/2 and j = 3/2 based bands with the 
same m = ± 1 / 2 . We have argued that the strong renor¬ 
malization of the SOC near the MIT as well as the full 
Hund’s coupling and pair hopping are essential for the 
emergence of the EXI state. 

Many experimental and theoretical studies support 
that the 5 d transition metal Ir oxides such as Sr 2 lr 04 and 
Ba 2 lr 04 are the j = 1/2 AF insulators [14-24]. These 
materials are characterized with the large SOC and mod¬ 
erate Coulomb interactions, thus in good qualitative ac¬ 
cordance with our phase diagram. Possible materials for 
the excitonic insulator found here should be characterized 
with i) larger Coulomb interactions, ii) relatively small 
SOC, and iii) the low spin configuration with no e g or¬ 
bitals involved. For 4 d electron systems, relatively large 
Coulomb interactions and relatively small SOC are ex¬ 
pected, while 3d electron systems exhibit large Coulomb 
interactions and very small SOC with the high or in¬ 
termediate spin configuration. We expect that the ideal 
systems for the stable excitonic insulator would be some¬ 
where between 3d and 4 d electron systems. This can 
be achieved in experiments by, for example, 3d and 4 d 
transition metal intermixed oxides or 4 d transition metal 
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oxide thin films grown on a compressive substrate with 
large distortion of oxygen octahedra. Further experimen¬ 
tal studies in this direction is highly desired. 
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